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We investigate the lepton flavor violating (LFV) rare decays in the supersymmetric minimal seesaw 
model in which the Frampton-Glashow-Yanagida ansatz is incorporated. The branching ratio of fi — ► 
cy is calculated in terms of the Snowmass Points and Slopes (SPS). We find that the inverted mass 
hierarchy is disfavored by all SPS points. In addition, once the ratio of BR(r — » fx'y) to BR(/i — » ej) 
is measured, one may distinguish the normal and inverted neutrino mass hierarchies, and confirm 
the masses of heavy right-handed Majorana neutrinos by means of the LFV processes and the 
thermal leptogenesis mechanism. It is worthwhile to stress that this conclusion is independent of 
the supersymmetric parameters. 
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O '■ I. INTRODUCTION 

<D ■ 

Recent solar pj, atmospheric Q, reactor Q and accelerator Q neutrino oscillation experiments have provided us with 
ON ', very robust evidence that neutrinos are massive and lepton flavors are mixed. The canonical seesaw mechanism Q 
• gives a very simple and appealing explanation of the smallness of left-handed neutrino masses - it is attributed to 
the largeness of right-handed neutrino masses. The existence of neutrino oscillations implies the violation of lepton 
flavors. Hence the lepton flavor violating (LFV) decays in the charged-lepton sector, such as /i — > e + 7, should 
also take place. They are unobservable in the Standard Model (SM), because their decay amplitudes are expected 
, to be highly suppressed by the ratios of neutrino masses (m i < 1 eV) to the W-boson mass (M w w 80 GeV). In 
the supersymmetric extension of the SM, however, the branching ratios of such rare processes can be enormously 
enlarged. Current experimental bounds on the LFV decays /i^e + 7, T^e + 7 and r — > \i + 7 are@ 
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The sensitivities of a few planned experiments Q may reach BR(/i — > ej) < 1.3 x 10 13 , BR(r — > ej) ~ O(10 8 ) and 
BR(r^/i 7 ) -O(10~ 8 ). 

For simplicity, here we work in the framework of the minimal supergravity (mSUGRA) extended with two heavy 
right-handed Majorana neutrinos. Then all the soft breaking terms are diagonal at high energy scales, and the only 
source of lepton flavor violation in the charged-lepton sector is the radiative correction to the soft terms through the 
neutrino Yukawa couplings. In other words, the low-energy LFV processes L — ► ij + 7 are induced by the RGE effects 
of the slepton mixing. The branching ratios of L — > l i + 7 are given by [1, Q 
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where m and A Q denote the universal scalar soft mass and the trilinear term at A GUT , respectively. In addition, mg 
is a typical mass of superparticles, can be approximately written as (Io| 

mf w 0.5m 2 M 2 /2 (m 2 + 0.6M 2 /2 ) 2 (3) 

with M t , 2 being the gaugino mass; and 

Ca = Y.(MMM^) jk \n±^L (4) 
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with Aqut = 2.0 x 10 16 GeV to be fixed in our calculations. M D and M i (for i = 1, 2) represent the Dirac neutrino 
mass matrix and the masses of the heavy right-handed Majorana neutrinos, respectively. 

To calculate the branching ratios of 2 • — > l i + 7, we need to know the following parameters in the framework of the 
mSUGRA: M x > 2 , m , A , tan/3 and sign(/i). These parameters can be constrained from cosmology (by demanding 
that the proper supersymmetric particles should give rise to an acceptable dark matter density) and low-energy 
measurements (such as the process b — > s + 7 and the anomalous magnetic moment of muon — 2). Here we adopt 
the Snowmass Points and Slopes (SPS) listed in Table. 1. These points and slopes are a set of benchmark 
points and parameter lines in the mSUGRA parameter space corresponding to different scenarios in the search for 
supcrsymmetry at present and future experiments. 

Since the canonical seesaw models are usually pestered with too many parameters, specific assumptions have to 
be made for M D or M i so as to calculate the LFV rare decays 112 . In addition, the LFV processes have also been 
discussed in the supersymmetric minimal seesaw model (MSM) [3/ B EE d , in which only two heavy right-handed 
Majorana neutrinos are introduced EH ■ I n the MSM, all model parameters can in principle be fixed by use of the 
LFV rare decays and eletric dipole moment of the electron . Ibarra and Ross discuss that the LFV processes may 
constrain the masses of right-handed heavy majorana neutrinos [TJ]. In the Frampton-Glashow-Yanagida [TtJ (FGY) 
ansatz, Raidal and Strumia have given analytic approximations to the LFV rare decays for the normal hierarchy 
case [l5|. As shown in the literature the FGY ansatz only includes three unknown parameters: M l5 M 2 and the 
smallest mixing angle Z . In this paper, we are going to numerically compute the LFV processes in terms of 9 Z for both 
normal and inverted hierarchies in the FGY ansatz. In the following parts, we shall show that the branching ratio of 
/i — » e7 may distinguish the neutrino mass hierarchies [42l|. In addition, once the ratio of BR(t — > fij) to BR(/i — > ey) 
is measured, one may distinguish the normal and inverted neutrino mass hierarchies, and confirm the masses of two 
heavy right-handed Majorana neutrinos by means of the LFV processes and the thermal leptogenesis mechanism pol]. 
It is worthwhile to stress that this conclusion is independent of the mSUGRA parameters. The remaining parts of this 
paper arc organized as follows. In Section II, we briefly describe the main features of the FGY ansatz and the thermal 
leptogenesis in the minimal supersymmetric standard model (MSSM). In Section III, the branching ratio of fi — > 
is computed in terms of the SPS. In Section IV, we numerically calculated the ratio of BR(r — > /17) to BR(/i — > e~f). 
Finally the summary are given in Section V. 



II. THE FGY ANSATZ AND THE THERMAL LEPTOGENESIS 



In the supersymmetric extension of the MSM, two heavy right-handed Majorana neutrinos (for i = 1,2) are 
introduced as the SU(2) L singlets. The Lagrangian relevant for lepton masses can be written as [21| 

- Acpton = hJi E R H i + hXu N K H 2 + I Af M R N R + h.c. , (5) 

where denotes the left-handed lepton doublet, while E R and N R stand respectively for the right-handed charged- 
lcpton and neutrino singlets. H 1 and H 2 (with hypcrcharges ±1/2) are the MSSM Higgs doublet superfields. After 
the spontaneous gauge symmetry breaking, one obtains the charged-lepton mass matrix M l = ViY, and the Dirac 
neutrino mass matrix M D = v 2 Y v . Here v i being the vev of the Higgs doublet H i (for i — 1,2). An important 
parameter (3 is defined by tan/3 = v 2 /v 1 or sin/3 = vijv with v ~ 174 GeV. The heavy right-handed Majorana 
neutrino mass matrix M R is a 2 x 2 symmetric matrix and M D is a 3 x 2 matrix. Without loss of generality, we 
work in the flavor basis where M ; and M R are both diagonal, real and positive; i.e., M x = Diag{m e , m„, m T } and 
M R = Diag{Af 1; M 2 }. The seesaw relation [| 

M v = -M D M^MS (6) 

remains valid. Note that this canonical seesaw relation holds up to the accuracy of 0(Mp/M R ) 22]. Since M R is 
of rank 2, M v is also a rank-2 matrix with \Det(M v )\ ~ m 1 m 2 m 3 = 0, where m i (for i — 1,2,3) are the masses 
of three light neutrinos. As for three neutrino masses m i , the solar neutrino oscillation data have set m 2 > m 1 [l|. 
Now that the lightest neutrino in the MSM must be massless, we are then left with either m l = (normal mass 
hierarchy) or m 3 = (inverted mass hierarchy). With the best-fit values A?7jg Un = m 2 — m\ = 8.0 x 10~ 5 eV 2 and 
Am 2 tm = \m\ — m 2 \ = 2.5 x 10 -3 eV 2 [23|, one can numerically calculate the neutrino masses [24| . 
In the FGY ansatz [ijj , M D is taken to be of the form 



M D = a 2 b 2 . (7) 
V bj 
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With the help of Eqs.(6) and (7), one may straightforwardly arrive at 
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Without loss of generality, one can always redefine the phases of left-handed lepton fields to make a l , b 2 and b 3 real 
and positive. In this basis, only a 2 is complex and its phase (f> = arg(a 2 ) is the sole source of CP violation in the model 
under consideration. Because a x , b 2 and b 3 of M D have been taken to be real and positive, M v may be diagonalized 
in a more general way 
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where P l = iDiag{e M , e 1 ^ , e* 7 } and V is lepton flavor mixing matrix (25j parameterized as 
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with s x = sin9 x , c x = cos6 x and so on. 6 a 

view of the current experimental data, we have 9 X = anu u y = <±o i^uesi-ni vaiuesj as wen as u z <^ ru ai me »»7o 
confidence level (23[. It is worth remarking that there is only a single nontrivial Majorana CP-violating phase (a) in 
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for the normal mass hierarchy (mi = 0). Similar results can also be obtained for the inverted mass hierarchy (m 3 = 0) 
p^ |. but we do not elaborate on them here. Because | cos<5| < 1 must hold, we find 0.077 < s z < 0.086 (mi = 0) and 
0.0075 < s z < 0.174 (m 3 = 0). Hence, a measurement of the unknown angle 9 Z becomes crucial to test the model. 

In the flavor basis where the mass matrices of charged leptons and right-handed neutrinos are diagonal, one can 
calculate the CP asymmetry in the decays of the lighter right-handed neutrino and obtain [27| 



M 1 |(MJ 12 | 2 |(MJ 23 | 2 sin20 
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Then e 1 can result in a net lepton number asymmetry Y L . The lepton number asymmetr y Y, is eventually converted 
into a net baryon number asymmetry Y B via the nonperturbative sphaleron processes [28l l29f. which is given by 

Y B = — ^ = -cY L = -c- £l , 14 

where c = 28/79 « 0.35, = 228.75 is an effective number characterizing the relativistic degrees of freedom which 
contribute to the entropy s of the early Universe, and k accounts for the dilution effects induced by the lepton-number- 
violating wash-out processes. The dilution factor k can be figured out by solving the full Boltzmann equations (27j. 
however, we take the following approximate formula (30j: 
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with m 1 = (A/qM d ) 11 /M 1 . It is clear that e 1 and Y B only involve two unknown parameters: M 1 and #z.[43| A 
generous range 8.5 x 10~ n < Yb 5. 9. 4 x 10~ n has been drawn from the recent Wilkinson Microwave Anisotropy 
Probe (WMAP) observational data [II]. It has been shown that M x > 3.4 x 10 10 GeV (mi = 0) and M 1 > 2.5 x 10 13 
GeV (m 3 = 0) are required by the current observational data of Y B [2l[. Flavor effects in the mechanism of thermal 
leptogenesis have recently attracted a lot of attention[33l|. Because the r Yukawa coupling are in thermal equilibrium 
for 10 9 GeV < Mi < 10 12 GeV, the flavor issue should be taken into account in our model for the mi = case. 
Using the analytic approximate formulas in the Ref. [HI] , one may calculate the CP-violating asymmetry and the 
corresponding wash-out parameter for the t lepton doublet in the final states of N l decays. We find e^. 1 = and 
= because of MD31 = 0, namely the N\ decays involving the t lepton doublet don't contribute to the baryon 
asymmetry Yb . On the other hand, the CP-violating asymmetry , which is produced by the out-of-equilibrium 
decays of N 2 , vanishes. Therefore, such flavor effects may be negligible in our paper. 

Note that there is in general a potential conflict between achieving successful thermal leptogenesis and avoiding 
overproduction of gravitinos in the MSM with supersymmetry [3o| unless the gravitinos are heavier than ~ 10 TeV 
[3^ ]. If the mass scale of gravitinos is of 0(1) TeV, one must have Mi < 10 s GeV. This limit is completely disfavored 
in the FGY ansatz. Such a problem could be circumvented in other supersymmetric breaking mediation scenarios 
(e.g., gauge mediation [lijl) or in a class of supersymmetric axion models [13], where the gravitino mass can be much 
lighter in spite of the very high reheating temperature. For simplicity, we choose Y B = 9.0 x 10~ n as an input 
parameter in this paper. Hence, one may analyze the dependence of Mi on 9 Z from the successful leptogenesis. 

III. DISTINGUISH THE NEUTRINO MASS HIERARCHY 

With the help of Eqs. (4) and (7), | | 2 can explicitly be written as 



\C 12 \ 2 = K| 2 |a 2 | 2 (ln^ 

|2 



2 



ic 13 r = 0, 



A, 



l^ 23 r = l&af \h\ 2 J ■ ( 16 ) 

Because of |C 13 | 2 = 0, we are left with BR(r — > ey) = 0. If BR(r — ► &y) 7^ is established from the future 
experiments, it will be possible to exclude the FGY ansatz. Using Eq. (13), we reexpress Eq. (16) as 

|C 12 | 2 = Ml \{M v ) l2 \ 2 (m%^ 2 

|C 23 | 2 = Af 2 |(A/J 23 | 2 (ln%^) 2 • (17) 

As shown in Section II, M 1 may in principle be confirmed by leptogenesis for given values of sin# z , but M 2 is entirely 
unrestricted from the successful leptogenesis with M 2 3> Mi . 

We numerically calculate BR(/i — > ey) for different values of sin^ by using the SPS points [13] ■ The results are 
shown in Fig. 1. Since the SPS points la and lb (or Points 2 and 3) almost have the same consequence in our 
scenario, we only focus on Point la (or Point 3). When sin# z — > 0.077 or sin6 z — > 0.086, the future experiment 
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is likely to probe the branching ratio of \x — > e + 7 in the m 1 = case. The reason is that sin# 2 — > 0.077 (or 
sin0 z — > 0.086) implies — > — 7r/2 (or </> — > 0). Furthermore, the successful leptogenesis requires a very large Mj due 
to s 1 oc sin2</>. It is clear that all SPS points arc all unable to satisfy BR(^t — > e7) < 1.2 x 10~ n in the m 3 = 
case. Therefore, we can exclude the m 3 = case when the SPS points are taken as the mSUGRA parameters. When 
s'm0 z ~ 0.014, BR(/i — > 67) arrives at its minimal value in the m 3 = case. For the SPS slopes, larger M x , 2 yields 
smaller BR(/Lt — > e7). We plot the numerical dependence of BR(/i — > e7) on Mjy 2 in Fig. 2, where we have adopted 
the SPS slope 3 and taken 300 GeV < M 1/2 < 1000 GeV. We find that M 1/2 > 474 GeV (or M 1/2 > 556 GeV) can 
result in BR(/t -» e7) < 1.2 x lO" 11 for sin0 z = 0.014 (or sin0 z = 0.1). For all values of M 1/2 between 300 GeV 
and 1000 GeV, BR(/j — > 57) is larger than the sensitivity of some planned experiments, which ought to examine the 
m 3 = case when the SPS slope 3 is adopted. The same conclusion can be drawn for the SPS slopes la and 2. In 
view of the present experimental results on muon g fl — 2, one may get M 1 , 2 < 430 GeV for tan [3 = 10 and A = p[39|. 
implying that the m 3 = case should be disfavored. 



IV. DETERMINE THE HEAVY MAJORANA NEUTRINO MASSES 



In this paper, an important conclusion is the masses of two heavy right-handed Majorana neutrinos (Mi and M%) 
can be derived through the leptogenesis and the LFV rare decays. With the help of Eqs.(2) and (17), one can obtain 

_ BR(r -> fj/y) _ Mj |(M„) 23 | 2 [ln(A GUT /M 2 )] 2 
' BR0we 7 ) Ml |(M„) 12 | 2 [ln(A GUT /Mi)] 2 ' 

Since the successful leptogenesis can be used to fix M ly a measurement of the above ratio will allow us to determine 
or constrain M 2 . On the other hand, we may calculate the ratio R by inputting the appropriate Mi- It is worthwhile 
to remark that the ratio R is independent of the mSUGRA parameters. To illustration, we show the numerical results 
of R as functions of sm6 z in Fig. 3(a) for the mi = case and Fig. 3(b) for the m.3 = case, respectively. In the 
m 3 = case, when sin 6^ — ► 0.0075, the phase 0^0 which implies Mi and M 2 may be larger than the GUT scale 
A GUT = 2.0 x 10 16 GeV. So we demand that Mi and M 2 must be less than the GUT scale. Below the scale A GUT , 
Mf (ln(A GUT /M 2 )) 2 obtain its maximum value at 

M 2 = A GUT /e = 7.4 x 10 15 GeV . (19) 

Consequently, we can derive the maximum value of R with M 2 = 7.4 x 10 15 GeV. In Fig. 3(b), the numerical 
results at the left side of the red (green) dash-dot line are "nonphysical" , because these results are corresponding to 
M 2 > 7.4 x 10 15 GeV for the M 2 = 10Mi (M 2 = 50M X ) case. We can obtain the ratio M 2 /M x once R is measured, 
furthermore calculate the heavy neutrino mass M 2 . From Fig. 3, one may derive R < 2 x 10 9 for the mi = case 
and R < 8 x 10 3 for the 7713 = case, respectively. If the future experiments prove R > 8 x 10 3 , the inverted mass 
hierarchy case can be excluded. It is worthwhile to stress that the above conclusions are independent of the mSUGRA 
parameters. 

Finally, let us comment on the Mi 3> M 2 case [l5[. In this case, the CP-violating asymmetry e 2 , which is produced 
by the out-of-equilibrium decay of N 2 , may finally survive. In order to producing the positive cosmological baryon 
number asymmetry (Yb > 0), 5 of Eq.(ll) must take - sign. The flavor effects may be neglected since the successful 
thermal leptogenesis requires M 2 > 10 12 GeV (mi = 0) and M 2 > 4.4 x 10 13 GeV (m 3 = 0). For Mi > 5M 2 , we 
calculate the branching ratio of fx — > ey, find the inverted hierarchy case is disfavored for all SPS points and slopes. 
For the normal hierarchy case, on may distinguish the Mi M 2 case (R > 14) and the Mi 3> M 2 case (R < 24) 
in terms of the ratio R. In addition, the masses of heavy neutrinos may be quasi-degenerate [40| . In the Mi ~ M 2 
case, the so-called resonant leptogenesis may occur |4lj . Such a scenario could allow us to relax the lower bound on 
the lighter right-handed Majorana neutrino mass. Therefore, the branching ratios of /j, — > e7 and t — » yU7 may be 
much less than the sensitivities of a few planned experiments. R ~ |(M^) 23 | 2 /|(M !/ )i 2 | 2 can be derived from Eq.(18). 
Furthermore, one may obtain 0.04 < R^ 1 < 0.07 (normal hierarchy) and 0.00045 < R^ 1 < 0.22 (inverted hierarchy). 



V. SUMMARY 



We have analyzed the LFV rare decays in the supersymmetric version of the MSM in which the FGY ansatz is 
incorporated. In this scenario, there are only three unknown parameters 9 Z , Mi and M 2 . The successful leptogenesis 
may fix Mi for given values of 8 Z . We have numerically calculated the branching ratio of /i — > e^j for different values 
of sm9 z by using the SPS Points. Then one can find the inverted mass hierarchy case is disfavored by all SPS points. 
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For the SPS slopes, the planned experiments may measure /i — > ej in the iris = case. In addition, once the ratio of 
BR(t — * fij) to BR(/j —> &y) is measured, we may also distinguish the neutrino mass hierarchies and fix the masses 
of two heavy right-handed Majorana neutrinos by means of the LFV processes and the leptogenesis mechanism. It is 
worthwhile to stress that this conclusion is independent of the mSUGRA parameters. 
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Point 


M 1/2 


m 


A 


tan p 


Slope 


la 


250 


100 


-100 


10 


mo = — Aq — OAAlj/2, M-i/2 varies 


lb 


400 


200 





30 


2 


300 


1450 





10 


m = 2Mx/n + 850 GeV, M 1/2 varies 


3 


400 


90 





10 


m = 0.25M 1/2 - 10 GeV, M 1/2 varies 


4 


300 


400 





50 


5 


300 


150 


-1000 


5 





TABLE I: Some parameters for the Snowmass Points and Slopes (SPS) in the mSUGRA. The masses are given in GeV. fi 
appeared in the Higgs mass term has been taken as /i > for all SPS. 
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FIG. 1: Numerical illustration of the dependence of BR(/i — > e^y) on sin# z : (a) in the = case; and (b) in the m 3 = 
case. The black solid line and black dash-dot line denote the present experimental upper bound on and the future experimental 
sensitivity to BR(/j — > &y), respectively. 
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FIG. 2: Numerical illustration of the dependence of BR(/i — > e*y) on M x / 2 f° r SPS slope 3 in the m 3 = case. The black 
solid line and black dash-dot line denote the present experimental upper bound on and the future experimental sensitivity to 
BR(/j — > e7), respectively. 
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FIG. 3: Numerical illustration of the dependence of BR(r — > /ry)/BR(/i — > ej) on sin# z : (a) in the mi = case; and (b) in 
the 7713 = case. 



